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ABSTRACT 

The  performance  of  three  recently  proposed  second-order  closure  models  is  tested  in 
benchmark  turbulent  shear  flows.  Both  homogeneous  shear  flow  and  the  log-layer  of  an 
equilibrium  turbulent  boundary  layer  are  considered  for  this  purpose.  An  objective  analysis 
of  the  results  leads  to  an  assessment  of  these  models  that  stands  in  contrast  to  that  recently 
published  by  other  authors.  A  variety  of  pitfalls  in  the  formulation  and  testing  of  second- 
order  closure  models  are  uncovered  by  this  analysis. 
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No.  NAS1-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  Virginia  23681-0001. 


INTRODUCTION 


The  need  for  advanced  turbulence  models  to  reliably  compute  the  complex  aerodynamic 
flows  of  technological  interest  has  led  to  a  resurgence  of  interest  in  second-order  closure  mod¬ 
els.  Consequently,  the  recent  papers  by  Shih,  Chen  and  Lumley1  and  Shih  and  Lumley2, 
which  reported  tests  of  second-order  closure  models  in  turbulent  shear  flows,  attracted  our 
attention.  In  these  papers,  results  were  presented  which  appear  to  indicate  that  the  Shih- 
Lumley  model3  performs  better  than  other  recently  proposed  second-order  closures  in  homo¬ 
geneous  shear  flow  as  well  as  in  other  boundary-free  turbulent  shear  flows.  However,  our  own 
comparative  studies  of  second-order  closure  models  in  benchmark  turbulent  shear  flows  have 
yielded  a  different  picture.  The  purpose  of  the  current  paper  is  to  present  these  alternative 
results  for  comparison. 

The  predictions  of  three  second-order  closure  models  recently  proposed  by  Shih  and 
Lumley2,3,  Fu,  Launder  and  Tselepidakis4  and  Speziale,  Sarkar  and  Gatski5  will  be  com¬ 
pared  in  two  benchmark  turbulent  flows:  homogeneous  shear  flow  and  the  log-layer  of  an 
equilibrium  turbulent  boundary  layer.  These  flows  are  selected  since  the  former  constitutes 
a  basic  building-block  free  turbulent  shear  flow  whereas  the  latter  serves  as  a  cornerstone  for 
the  calculation  of  practical  wall- bounded  turbulent  flows  of  engineering  interest.  Particular 
attention  will  be  paid  to  evaluating  the  ability  of  each  model  to  accurately  predict  the  equi¬ 
librium  values  for  the  Reynolds  stress  anisotropies.  However,  for  the  case  of  homogeneous 
shear  flow,  model  predictions  for  time  evolving  Adds  will  also  be  compared.  Objective  means 
for  evaluating  the  performance  of  the  models  will  be  provided  and  pitfalls  in  the  formulation 
and  evaluation  of  models  are  uncovered  that  have  led  to  previously  published  assessments 
that  are  misleading. 

THE  TURBULENT  SHEAR  FLOWS  TO  BE  CONSIDERED 

We  will  consider  incompressible  turbulent  shear  flows  with  the  mean  velocity  gradient 
tensor 

||  =  SM*  (i) 

where  Sij  is  the  Kronecker  delta  and  S  is  shear  rate.  In  homogeneous  shear  flow,  the  shear  rate 
S  is  constant  and  is  applied  in  an  unbounded  flow  domain  yielding  spatially  homogeneous 
turbulence  statistics.  For  this,  as  well  as  any  homogeneous  turbulence,  the  Reynolds  stress 
tensor  r,j  is  a  solution  of  the  transport  equation6 
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where  t,3  =  u.tij  and 
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are,  respectively,  the  pressure-strain  correlation  and  dissipation  rate  tensor  (here,  p  is  the 
fluctuating  pressure,  ti,  is  the  fluctuating  velocity,  u  is  the  kinematic  viscosity,  and  an  overbar 
represents  an  ensemble  mean).  Since 


(4) 


where  II, ■  j  =  —  o£ij  given  that  o^ij  is  the  deviatoric  part  of  the  dissipation  rate  tensor 

(e  =  closure  is  achieved  once  models  for  n,j  and  s  are  provided.  In  most  existing 

second-order  closure  models,  is  neglected  while  II  tJ-  and  e  are  modeled  in  the  general 
form: 


n,;  =  eAi(6)  +  KMtjki(b)~ 

e  &Ui  e2 
e  =  —Ce\—Tu— - Crt— 


KxldXi  elK 


(5) 

(6) 


where 


,  Ta  ~  I KStj 

b,j  ~  2  K 


(7) 


are  the  turbulent  kinetic  energy  and  Reynolds  stress  anisotropy  tensor,  respectively.  Here, 
Aij  and  Miju  are  dimensionless  tensor  functions  of  and  possibly  the  turbulence  Reynolds 
number  Rt  =  K2 / i/e;  Cel  and  CC2  are  either  constants  or  functions  of  the  second  and  third 
invariants  (//,///)  of  6,7  as  well  as  Rt.  The  full  form  of  the  three  models  to  be  considered 
-  the  Shih-Lumley  (SL)  model,  the  Fu,  Launder  and  Tselepidakis  (FLT)  model  and  the 
Speziale,  Sarkar  and  Gatski  (SSG)  model  -  are  provided  in  the  Appendix. 

For  homogeneous  shear  flow,  each  of  the  models  -  consistent  with  physical  and  numerical 
experiments  -  predict  that  the  anisotropy  tensor  6,,  and  shear  parameter  SK/e  achieve  equi¬ 
librium  values  that  are  independent  of  the  initial  conditions  (see  Tavoularis  and  Corrsin7, 
Rogers,  Moin  and  Reynolds8  and  Tavoularis  and  Karnik9).  This  equilibrium  state  is  associ¬ 
ated  with  solutions  where  =  0  or  equivalently, 
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where  V  =  —Tijdvi/dxj  is  the  turbulence  production.  The  substitution  of  (1)  and  (8)  into 
(2)  yields  the  system  of  algebraic  equations  (see  Abid  and  Speziale10) 
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where  Il’y  =  II ij/SK  is,  for  the  turbulent  shear  flows  to  be  considered,  a  function  of  Tij/K 
and  V ft  whose  specific  form  depends  on  the  pressure-strain  model  chosen.  In  deriving  (9), 
the  identity 

?-£(“)  w> 

has  been  used.  Once  the  equibbrium  value  of  "P/e  is  specified  (and  Eq.  (7)  is  utilized),  it 
is  straightforward  to  obtain  the  equilibrium  values  of  b,3  from  a  numerical  solution  of  (9). 
These  equilibrium  values  are  determined  exclusively  by  the  pressure-strain  model. 

In  the  logarithmic  region  of  an  equilibrium  turbulent  boundary  layer,  the  mean  velocity 
gradient  tensor  is  of  the  general  form  (1)  and  there  is  a  production-equals-dissipation  equi¬ 
librium  where  the  turbulent  diffusion  terms  vanish  in  the  Reynolds  stress  transport  equation. 
Consequently,  (9)  yields  the  equilibrium  Reynolds  stress  anisotropies  for  the  log-layer  when 
P/e  =  1  (see  Abid  and  Speziale10).  In  contrast  to  this,  homogeneous  shear  flow  achieves 
an  equilibrium  state  where  P/e  >  1  (physical  and  numerical  experiments7-9  have  indicated 
that  1.4  <  P/e  <  1.8). 


DISCUSSION  OF  RESULTS 

A  comparison  of  the  results  predicted  by  the  pressure- strain  models  of  Shih  and  Lumley2,3, 
Fu,  Launder  and  Tselepidakis4  and  Speziale,  Sarkar  and  Gatski5  will  now  be  made.  In 
Figures  1(a)  and  1(b),  the  model  predictions  for  the  norm  of  the  slow  and  rapid  parts  of  the 
pressure-strain  correlation  are  compared  with  the  DNS  results  of  Rogers,  Moin  and  Reynolds8 
for  homogeneous  shear  flow  (II, ^  =  Ilj^  +  Iljj^  where  the  superscripts  (5)  and  (R)  denote, 
respectively,  the  slow  and  rapid  parts;  (ILjlL,)1/2  denotes  the  L2  norm  of  II, j)).  These 
results  are  similar  to  those  used  by  Shih  and  Lumley2  to  conclude  that  the  Shih-Lumley 
model  performs  the  best  of  these  three  models  and  that  the  SSG  model  performs  poorly. 
Such  a  conclusion  is  highly  misleading.  The  SSG  model  is  a  model  for  the  total  pressure- 
strain  correlation  and  not  for  its  separate  slow  and  rapid  parts.  The  standard  hierarchy  of 
pressure-strain  models  (5)  are  only  theoretically  justified  for  homogeneous  turbulent  flows 
that  are  near  equilibrium  (as  shown  by  Speziale,  Gatski  and  Sarkar11 ,  both  the  slow  and  rapid 
parts  of  the  pressure-strain  correlation  depend  nonlinearly  on  the  mean  velocity  gradients  at 
retarded  times  in  non-equilibrium  turbulent  flows).  When  (5)  is  thought  of  as  the  simplified 
equilibrium  form  of  a  more  general  pressure-strain  model  it  becomes  ambiguous  as  to  which 
part  of  (5)  represents  the  slow  pressure-strain  and  which  part  of  (5)  represents  the  rapid 
pressure-strain.  It  can  only  be  said  definitively  that  e Aij  is  the  slow  pressure-strain  in  the 
limit  of  relaxational  turbulent  flows  and  that  KMijkidvk/dxi  is  the  rapid  pressure-strain  in 
the  rapid  distortion  limit.11  This  ambiguity  causes  no  problem  since  only  the  total  pressure- 
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strain  correlation  is  needed  in  (2)  for  the  calculation  of  the  Reynolds  stresses.  When  the 
model  predictions  for  the  norm  of  the  total  pressure-strain  correlation  are  compared  with 
the  DNS  of  Rogers  et  al.s  for  homogeneous  shear  flow  a  rather  different  picture  emerges  as 
can  be  seen  in  Figure  1(c).  The  SSG  model  performs  as  well,  if  not  better,  than  the  SL  and 
FLT  pressure-strain  models.  Furthermore,  it  was  shown  by  Speziale,  Sarkar  and  Gatski11 
that  none  of  the  other  models  are  capable  of  predicting  the  individual  slow  or  rapid  parts  of 
the  pressure-strain  correlation  for  a  wide  range  of  homogeneous  turbulent  flows.  The  results 
presented  in  Shih  and  Lumley2  are  misleading  in  this  regard. 

Comparisons  such  as  those  shown  in  Figure  1  are  not  very  helpful  for  determining  the 
predictive  capabilities  of  a  pressure-strain  model  in  turbulent  shear  flows.  A  better  test  for 
gauging  the  performance  of  a  model  is  to  determine  its  ability  to  predict  accurate  equilib¬ 
rium  values  for  the  Reynolds  stress  anisotropies.  These  are  the  crucial  physical  quantities 
in  homogeneous  shear  flow  that  are  independent  of  the  initial  conditions  and,  therefore, 
repeatable.  Over  thirty  independent  test  runs  of  homogeneous  shear  flow  in  reported  phys¬ 
ical  and  numerical  experiments7-9  have  yielded  equilibrium  values  for  the  Reynolds  stress 
anisotropies  that  lie  within  10%  of  one  another.  On  the  other  hand,  the  time  evolutions  of 
the  individual  Reynolds  stresses  vary  by  factors  of  3  or  4  depending  on  the  initial  conditions 
of  the  test  case. 

In  Table  1,  the  equilibrium  Reynolds  stress  anisotropies  predicted  by  the  SL,  FLT  and 
SSG  models  are  compared  with  the  experimental  results  for  homogeneous  shear  flow  for 
P/e  =  1.5.  This  value  of  P/e  is  chosen  since  it  is  the  average  value  obtained  from  the 
most  recent  experiments.  The  SSG  model  is  in  closest  agreement  with  the  experimental 
data  whereas  the  SL  model  exhibits  the  largest  discrepancies.  The  same  trend  is  exhibited 
for  other  equilibrium  values  of  P/e  in  the  experimental  range  of  1.4  -  1.8  (see  Abid  and 
Speziale10  for  more  details). 

As  discussed  earlier,  the  equilibrium  Reynolds  stress  anisotropies  for  the  log-layer  of  a 
two-dimensional  turbulent  boundary  layer  can  be  obtained  from  (9),  after  setting  P/e  =  1. 
The  results  obtained  from  the  SL,  FLT  and  SSG  models  are  compared  with  experimental 
data12,13  in  Table  2.  The  results  obtained  are  similar  to  those  obtained  for  homogeneous 
shear  flow:  the  SSG  model  is  in  the  closest  agreement  with  the  experimental  data  while  the 
SL  model  exhibits  the  largest  deviations.  In  regard  to  the  latter  point,  it  is  surprising  how 
poorly  the  SL  model  performs  in  its  prediction  of  the  normal  Reynolds  stress  anisotropies. 

In  Figure  2,  the  model  predictions  for  the  time  evolution  of  the  turbulent  kinetic  energy 
in  homogeneous  shear  flow  are  shown  for  two  test  cases:  the  large-eddy  simulation  (LES) 
of  Bardina  et  al.14  and  the  DNS  of  Rogers  et  al.8  (run  C128X).  Here,  the  dimensionless 
turbulent  kinetic  energy  and  dimensionless  time  are  given  by  K*  =  K/Kq  and  t*  =  St, 
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respectively.  These  results  exhibit  the  same  trend  as  before:  the  SSG  model  performs  the 
best  and  the  SL  model  has  the  largest  deviations.  The  low  growth  rate  predicted  by  the  SL 
model  arises  from  its  underprediction  of  the  equilibrium  value  of  bi2  (see  Table  1). 

This  brings  us  to  the  basic  question  as  to  why  the  SL  model  performs  poorly  in  these 
simple  shear  flows.  As  discussed  in  Shih  et  al.1 ,  the  SL  model  was,  to  a  large  extent,  calibrated 
based  on  realizability  constraints.  It  was  recently  shown  by  Speziale,  Abid  and  Durbin15  that 
the  SL  model  yields  unrealizable  results  in  homogeneous  shear  flow  due  to  an  error  in  their 
analysis.  In  Figure  3,  the  time  evolution  of  the  invariant  function  F  =  1  +  97/  +  27/77 
predicted  by  each  model  in  homogeneous  shear  flow  is  shown  for  the  anisotropic  initial 
conditions:  bn  =  —0.32,  622  =  &33  =  0.16,  b12  =  623  =  b13  =  0,  and  SK/e  —  15.  For 
realizable  turbulence,  we  must  have  F  >  0.  It  is  clear  that  for  these  initial  conditions,  the 
FLT  and  SSG  models  yield  realizable  solutions  whereas  the  SL  model  yields  unrealizable 
results!  Hence  the  primary  theoretical  constraint  by  which  the  SL  model  was  formulated  is 
in  error. 

In  conclusion,  it  must  be  said  that  the  good  performance  of  a  model  in  these  basic 
turbulent  shear  flows  is  no  guarantee  that  it  will  perform  well  in  more  complex  turbulent 
flows.  However,  these  simple  test  flows  do  bear  directly  on  how  well  a  model  will  perform 
in  equilibrium  turbulent  boundary  layers  which  form  a  cornerstone  for  many  engineering 
applications.  A  model  that  cannot  predict  these  simple  test  cases  accurately  should  be 
abandoned  for  use  in  more  complex  engineering  flows. 
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APPENDIX 


The  detailed  form  of  the  turbulence  models  considered  in  the  paper  are  as  follows: 


Shih  &  Lumley  Model 


Djj  —  —Ciebij  + -KSij  +  12a5K  (bikSjk  +  bjtSik 


2  _  \  4 

—  gh klSkiSijJ  4-  -(2  —  7as)K(b{kUjk  +  fyfcu;,*) 


+-^K(bnbimSjm  +  bj[btmSim  —  2bikSkibij 


-3bkiSkibij)  +  -K(bublmUjm  +  bjiblmu>im) 


Ci=  2  +  ^  exp(—7 .77  /  yfRet){72/  y/Ret  +  80.1  ln[  1  +  62.4(-//  +  2.3///)]}  (A2) 


F  =1  +  911  +  27  III 

II  =  ~2  bjbij,  III  =  ^bijbjkbkt 
Xe,  = 


°‘=To(1  +  lFk) 

T 

C.1  =  1.20,  Ce2  =  -  +  0.49  exp(-2.83/ v^)[l  -  0.33Zn(l  -  55//)] 


/\t,  Launder  &  Tselepidakis  Model 


Du  —  —Ci ebij  +  C2e  (btkbkj  —  - 

+  +  l-2/f  (bikSjk  +  bjkSik  -  -buSkiSij'j 

26  4 

+Y^ZC(6jfcu7jk  -|-  bjk&ik)  +  -K(bikbkiSji 
+bjkbkiSu  -  2bikSkibij  —  SbkiSkibij) 

+  (bikbu&ji  +  bjkbki&n)  -  —K  [81 1 (bikUjjk 
+bjk&ik)  +  12  (bikbu<3imbmj  +  bjkbkivimbmi)] 


12.  Klebanoff,  P.  S.,  “Characteristics  of  Turbulence  in  a  Boundary  Layer  with  Zero  Pressure 
Gradient,”  NACA  TN  5175, 1956. 

13.  Laufer,  J.,  “Investigation  of  Turbulent  Flow  in  a  Two-Dimensional  Channel,”  NACA 
Technical  Report  1058,  1951. 

14.  Bardina,  J.,  Ferziger,  J.  H.  and  Reynolds,  W.  C.,  “Improved  Turbulence  Models  Based 
on  Large-Eddy  Simulation  of  Homogeneous  Incompressible  Turbulent  Flows,”  Stanford 
University  Technical  Report  TF-19,  1983. 

15.  Speziale,  C.  G.,  Abid,  R.  and  Durbin,  P.  A.,  “New  Results  on  the  Realizability  of 
Reynolds  Stress  Turbulence  Closures,”  ICASE  Report  98-76,  NASA  Langley  Research 
Center,  1993. 
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Equilibrium 

Values 

SL  Model 

FLT  Model 

SSG  Model 

Experimental 

Data 

bn 

0.105 

0.177 

0.214 

0.20  to  0.21 

b\i 

-0.121 

-0.157 

-0.163 

-0.14  to  -0.16 

&22 

-0.107 

-0.122 

-0.140 

-0.14  to  -0.15 

633 

0.002 

-0.055 

-0.074 

-0.05  to  -0.07 

Table  1.  Equilibrium  Reynolds  stress  anisotropies  in  homogeneous  shear  flow:  Comparison 
of  the  model  predictions  with  physical  and  numerical  experiments.7-9 
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Equilibrium 

Values 

SL  Model 

FLT  Model 

SSG  Model 

Experimental 

Data 

611 

■ 

0.141 

0.201 

0.20  to  0.23 

b\2 

-0.116 

-0.162 

-0.160 

-0.14  to  -0.16 

-0.082 

-0.099 

-0.126 

-0.13  to  -0.15 

&33 

0.003 

-0.042 

-0.075 

-0.05  to  -0.10 

Table  2.  Reynolds  stress  anisotropies  in  the  log-layer  of  an  equilibrium  turbulent  boundary 
layer:  Comparison  of  the  model  predictions  with  experiments.12-13 
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Figure  1.  Comparison  of  the  model  predictions  for  the  norm  of  the  pressure- strain  correlation 
with  the  direct  numerical  simulation  of  Rogers  st  al.B  for  homogeneous  shear  flow  (run 
C128U):  (a)  slow  part,  (b)  rapid  part  and  (c)  total. 
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